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. We present an elementary model of the collapses and revivals in the visibility of the interference 

between two atomic Bose-Einstein condensates. We obtain different predictions of the revival times 
whether we conserve or break atom number conservation from the outset. The validity of Bose- 
broken symmetry can be tested by observations of these collapses and revivals. 
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The recent experimental realization of a weakly interacting Bose-Einstein condensate in an alkalic gas jj]-|3| has 
stimulated considerable theoretical work on the properties of these condensates. One question which has received 
much attention concerns the phase of the condensate and how it is established [^-|TT|] . A conventional approach 
is to invoke Bose-broken symmetry arguments [12] and select an arbitrary phase from an initial state which has a 



random phase. (It is not possible to measure the absolute phase of a condensate so that the phases we talk about 
are the relative phases between two condensates.) A condensate formed in the ground state of a trap will be in a 
number state, though we may not have knowledge of what this number is. A number state may be considered as a 
"-^J I continuous superposition of coherent states all with different phases; spontaneous symmetry breaking then selects just 
one (arbitrary) phase. Recently it has been shown that a relative phase may be established between two condensates 
initially in number states by measurements This leads to the conclusion that considering the condensate to 

be in either a coherent state or a number state are equivalent. Note that whereas atom number is conserved in the 
number state case, an initial coherent state is a Poissonian superposition of number states and hence the atom number 
conservation is broken. 

In this paper we consider the evolution of the visibility in the interference between two Bose-Einstein condensates 
in the presence of collisions. Collisions give rise to the collapse and revival of the macroscopic wave-function of small 
atomic condensates composed of 10 3 — 10 5 atoms |l3| , which results in collapse and revival in the visibility when two 
t— I ■ condensates are interfered. The central result of this paper is that the revival time is strongly dependent (a factor of 
2) on the initial state chosen for the condensate. This gives us a way to determine unambiguously the true quantum 
state of the condensate. We shall illustrate our results for three different initial states of the two condensates. In 
the first example we invoke Bose-broken symmetry and consider the product of two coherent states so that number 

■ conservation is broken. In the other two examples we consider an entangled state of the two condensates with fixed 
total number N so that no assumption of Bose-broken symmetry is invoked. The evolution of the condensates under 
the influence of collisions is studied for each of the above initial conditions and the collapse and revival times are 

' O I compared. 

Our model comprises two Bose-Einstein condensates which both occupy the ground-state of their respective traps 
and are described by the atom annihilation (creation) operators a (aJ) and b (b'). Atoms are released from each 
trap with momenta kiand \s.2 respectively, producing an interference pattern which enables a relative phase to be 
measured. 

■ The intensity of the atomic field is given by 

I(x,t)=I ( a+ (t) e* klOE + &+ (t) e <ka0t a (t) e^ 1 '* + b (i) e" ik2X ), (1) 



where Iq is the single atom intensity. Atoms within each condensate collide and this may be described by the 
Hamiltonian 



H = -hv 

2 A 



(2) 



where x 1S the collision rate between the atoms within each condensate. Cross-collisions between the two condensates, 
described by the term a'ab'b, are not included since they are dependent on the actual geometry of the physical 
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situation. The coefficient of this term could be anywhere between zero and h\ depending on the overlap between the 
two condensates. We consider the case where this overlap is small. 

Including the time dependence of a and b due to the collisions described by Eq. (||), we get for the intensity 

I (x,i) = I {(a f a) + {tfb\ + Itf exp i (a f a - tfb\ X t &) e" l0(x) + h.c.\ , (3) 

where tfi (x) = (k.2 — ki) • x. The third term in the expression for I (x,i) gives rise to interference fringes. However the 
interference pattern will be modified in time due to the decohering effects of the collisions. 

We shall now proceed to demonstrate how the interference pattern evolves for different initial states of the conden- 
sates. We will consider two classes of initial states, one which uses Bose-broken symmetry arguments and the other 
which does not. This Bose symmetry is associated with atom number conservation, so breaking this symmetry allows 
us to write the wave-function as a superposition of number states. This symmetry is deemed to be broken in the 
formation of the condensate. The uncertainty in the number leads to the wave-function possessing a definite phase 
via the number-phase uncertainty relation AnA0 ~ 1. Following this line of argument one can imagine an initial 
fixed number of atoms before the condensate is formed. This has indeterminate phase since it can be expressed as a 
continuous superposition of coherent states each possessing a particular phase. When the condensate is formed, the 
symmetry is broken with one of the coherent state selected, its phase becomes the phase of the condensate. In the 
spirit of Bose-broken symmetry we then take each condensate to be initially in a coherent state 

Wb) = \a) |/3) . (4) 

This yields for the intensity in Eq. (|J) 

I(x,t) = Iq \a\ 2 |l + exp 2|a| 2 (cosxi-l) cos [0 (x) - 0] | (5) 

where we have set the amplitude of the two condensates to be equal to maximize the visibility of the interference 
pattern. The relative phase between the two condensates is defined to be </>, so that the relationship between the 
complex amplitudes is (3 — ae" 1 *. The exponential term describes the time dependence of the visibility of the 
interference pattern. Inside this exponential, we have a periodic function of period 2tt/x, corresponding to revival 
times where the visibility is 1. This visibility suffers a minimum half-way between these revivals with a value of 
exp ^—4 |a| 2 ^ ; it varies smoothly in time from these minima to their local maxima. 

In our second class of initial states we will not use Bose-broken symmetry, so we conserve the total atom number 
N of the two condensates. We shall consider two entangled states of the condensates with fixed total number N. In 
the first example we consider the product of two coherent states \a) €5 \(3) projected onto a number state basis [ pr| . 
This basis is truncated to size N with equal amplitudes \a\ = \/3\ = y/ N/2. We can define a relative phase between 
the condensates by superposing number difference states. This entangled state is 
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\tp N ) = |fc) ® \N - k), (6) 

y/k\ (N-ky. 

where <j> is the relative phase between the condensates. Note how each entangled number state has fixed total atom 
number N and number difference 2k. The intensity for this initial state is 

J(x,t) =7oY{l + cos Ar - 1 X tcos[0(x)-0]}. (7) 

The visibility of the interference pattern is Jcos^^ 1 The parity of the total atom number N, whether it is odd or 
even, plays a role in the revivals. When N is odd, the cosine term is raised to an even power giving a revival period of 
7r/x since it is never negative. When N is even, the cosine term is raised to a odd power also giving a revival period 
of ir/x but with each alternate revival occurring with the phase shifted by tt radians. In both cases, the revival period 
is one half of the period predicted in the previous case where we used Bose-broken symmetry whatever the parity of 
N is. 

The significant difference between the two cases is that the total number N is not fixed for the case of Bose-broken 
symmetry whereas it is for the other case. The factor of two difference in the period can be explained by looking at 
the exponential term in Eq. (|^). Inside this exponential we have the atom number difference operator a^a— b^b which 
is quantized in units of 2 when the total number is fixed and units of 1 when it isn't fixed. Thus we have either an 
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exp (i2nxt) term or an exp(inxt) term, where n is an integer. This gives rise to the factor of two difference in the 
period. 

As a third example we shall consider an initial state which may be formed by quantum measurement. The state 
formed from two condensates initially in number states N/2 after m atoms have been detected is of the general form 
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Cfe|n — m + k, n — k), (8) 
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which is normalized so that the coefficients satisfy the relation ^ |cfc| 2 = 1, with the values of the c k depending on the 
actual sequence of measurements. Note that this prepared state has a fixed total number of atoms (2n — m atoms). 
In order to study the collapse and revival of this coherence we need to include the decohering effects of collisions. 
A Monte-Carlo wave-function method has been used to simulate numerically the time evolution between subsequent 
atom position measurements including the effects of collisions [|]|]]. The model we have used so far already include 
collisions so all that we need to do is to consider the prepared state \ f m ) as our third description of the wave- function 
for the two condensates. Using the expression for the intensity given by Eq. (||), we obtain 

/ (x,i) = Iq < 2n — m + c* k Ck-i y/ (n — k) (n — m + k — 1) exp [i (2k — m — 1) %i — i<fi (x)] + h.c. > . (9) 

I fc=l J 

Any collisional effects occurring during the state preparation, which we have not included here, will affect the coef- 
ficients of \iprn}- For a good state preparation in the sense that the entangled state possesses high coherence we will 
expect that the relative phases between the neighboring coefficients c& are very similar and in fact give a good estimate 
of the relative phase between the two condensates. Looking back at the coefficients previously defined in the second 
description of the wave- function, Eq. (^|), we see that the phase between the neighboring coefficients is identical for 
\<Pn)- By comparing the coefficients of the prepared state \(f m ) in numerical simulations where small collisional rates 
have been included (see reference ||), we know that there is a fairly well defined phase between neighboring entangled 
number states in \ip m ) so that we may write, to a good approximation, c^Ck-i = Ake 1 ^, where <f> is the relative phase 
between the condensates. The interference pattern then becomes 



I (x,<) = Jo \ 2ri — m + A k cos I( 2fc -m-l)xt + (l>-<f> (x)] \ , (10) 

where 



k=i 



A k = A ky /(n-k) (n-m + k-1), (11) 

which consists of a summation over cosines weighted by Ak with differing multiple frequencies of \t. Separating the 
time dependence from the phase 4> and the position <f> (x) terms in the cosine, we obtain 



J(x,t) 



2n — m + Ak cos [(2k — m — 1) %i] cos [<p (x) 



k=l 



(12) 



2k — m — 1 is odd (even) whenever m is even (odd). Since the initial total number is 2n which by definition is even, 
the parity of m determines the parity of the total atom number of the entangled state. When the parity is even, the 
intensity consists of a constant term plus a weighted sum over cosines with frequencies which are in multiples of 2\t 
and thus the period is tt/x- On the other hand, if the parity is odd we have a weighing over cosines with frequencies 
which are odd multiples of The "visibility" of the interference pattern I (t) would then be —1 or 1 for odd or 
even multiples of the revival time tt/x respectively. Visibility in the usual sense is always positive; the change in sign 
indicates a n phase shift in the interference pattern. Thus these two examples of entangled states which preserve 
atom number agree on the period and parity properties of the collapses and revivals. 

The visibilities of the interference patterns when Bose-broken symmetry is invoked, Eq. (^) and when it is not, 
Eq. (0) are shown in Fig. (|l|a) and (0b). We have not shown a graph of the third description because it is implicitly 
a measurement scheme which establishes phase and would give a result very similar to Eq. (Qb). In order to verify 
experimentally the dependence of the visibility on time shown in Fig. (Q), the following procedure may be employed. 
Initially a sequence of measurements of the output field is made in order to establish a relative phase between the 
condensates. The position of the peaks of this initial interference pattern is recorded. The condensates are then left 
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to evolve freely under the influence of the Hamiltonian Eq. (£|) . After a time t of free evolution, a second sequence of 
measurements of the output field is made. This establishes a new interference pattern with peaks located at positions 
which may differ from those of the initial pattern. The calculated visibility refers to the distribution of the position 
of the peaks of the new interference pattern relative to those of the initial pattern. As such it can only be determined 
by several pairs of measurement sequences which alternatively establish and probe the phase, separated by the time 
interval t. 

The validity of imposing Bose-broken symmetry from the outset as a means of describing the quantum dynamics of 
Bose-condensed systems can be tested by measuring the visibility of the interference pattern between two condensates 
as a function of time. For small atomic samples composed of 10 3 — 10 5 atoms the visibility undergoes a series of 
collapses and revivals. The magnitude of the revival period of the visibility will test whether the symmetry is in fact 
broken in the condensation process or whether a quantum state involving a fixed number of atoms is more appropriate. 
Indirect measurements are also possible from recent light scattering proposals suggested by Imamoglu and Kennedy 
and Javanainen |l5| using independent condensates coupled to a common excited state. The recent production 
of two overlapping condensates of 87 Rb in different spin states has increased the experimental feasibility of these 
proposals.. 
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FIG. 1. The visibility as a function of time for (a) two independent coherent states which invokes Bose-broken symmetry 
and (b) entangled coherent states which does not. We have used a total atom number of N — 100 corresponding to a square 
amplitude |a| 2 = 50 for the two independent coherent states. 
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